Controlling transition probability from matter-wave soliton to chaos 
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For a Bose-Einstein condensate loaded into a weak traveling optical superlattice it is demonstrated 
that under a stochastic initial set and in a given parameter region the solitonic chaos appears with 
a certain probability. Effects of the lattice depths and wave vectors on the chaos probability are 
investigated analytically and numerically, and different chaotic regions associated with different 
chaos probabilities are found. The results suggest a feasible method for eliminating or strengthening 
chaos by modulating the moving superlattice experimentally. 
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Many phenomena observed in Bose-Einstein conden- 
sates (BECs) are well modelled by nonlinear Schrodinger 
Equation (NLSE), also known as Gross-Pitaevskii equa- 
tion (GPE), which includes many fantastic nonlinear ef- 
fects, such as chaos [10, soliton[3j,kfl, and so on. Chaotic 
soliton behaviors, which are of particular interest, have 
been studied theoretically in a NLSE with a periodic 
perturbation 0, @, @]. Lately, there are a few remark- 
able works on chaotic dynamics of soliton in BEC sys- 
tems, including the chaos and energy exchange^, the 
discrete soliton and chaotic dynamics in an arra y o f 
BECs[9], and the bright matter- wave soliton collision [10|. 
Recently, increasing interest is excited in a BEC held 
in an optical superlattices for the periodic [ll], [l2j and 
quasiperiodic[H, [13, HU cases which are in close analo- 
gies with the fields of supercrystals and quasicrvstals [l6| . 
The BECs interacting with a traveling lattice have also 
been successfully treated experimentally [13, EH and 
theoretically [H S M that shows many fantastic results 
such as lensing eff ect Il8ll, gap soliton generationpjj], spa- 
tiotemporal chaos |20l l21j . and so on. From the above 
analyzes we get a physical motivation, namely using a 
moving optical superlattices to study the transitions from 
soliton to chaos in BEC matter. 

It is well known that for stochastic initial and bound- 
ary conditions and fixed parameters, chaos does not al- 
ways appear in a chaotic system[20|, but with certain 
probability. Chaos probability may play a significant role 
for the control of chaos. Many works have focused on 
suppressing chaos which results in zero chaos probabil- 
ity. Whereas, in some realistic applications such as secure 
communication based on chaos |22|. higher chaoticity is 
desired [23], which calls for higher chaos probability. In 
this paper, we show that in a BEC system perturbed by a 
moving optical superlattice consists of two lattices of dif- 
ferent depths and wave vectors, the superlattice can sep- 
arate the chaotic region into several parts with different 
chaos probabilities. Furthermore, for a fixed first lattice 
the adjustment to the secondary lattice could turn the 



chaos probability to zero or higher one. The results sug- 
gest a feasible method for eliminating or strengthening 
chaos experimentally. Such a method can be extended 
to the controls of the spatial chaos with zero traveling 
velocity and the temporal chaos in other systems. 

Primary and secondary moving optical lattice compose 
the superlattice [24| as the form 

V op (C) = Vt cos 2 (fcC) + V 2 cos 2 ^ + 0), (1) 

where we refer to V\ cos 2 '(kQ as the primary lattice with 
Vi and k corresponding to its depth and wave vector, 
and V2 cos 2 (-fk( + (j>) as the secondary one with V2 be- 
ing its depth, 7 the ratio of the two laser wave vectors 
and <f> the phase difference. The spatiotemporal variable 
£ = x + viX contains the velocity of the traveling lat- 
tice vl — S I '(2k) with S being the frequency difference 
between the two counter-propagating laser beams pro- 
ducing the first lattice. The laser frequencies and am- 
plitudes can be controlled independently by using the 
acousto-optic modulators [17j. 

Chaotic Solitons For the quasi-one-dimcnsional (ID) 
BEC held in superlattice potential V op (() the transverse 
wave function is approximately in ground state of a har- 
monic oscillator of frequency uj r , and the governing time- 
dependent GPE reads 
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dt 2m dx 2 
Here "J and m are the macroscopic quantum wave func- 
tion and mass, g\d = 2fjjjj r a s characterizes the ID inter- 
atomic interaction strength [25j with a s being the s-wave 
scattering length. In order to get a simple description, 
we investigate the traveling wave solution in the form[2(i| 
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where -R(C) and #(£) are real functions of £, a and j3 rep- 
resent two undetermined real constants. Inserting Eq.Q 
into Eq.([2]), we easily obtain two coupled equations 

d 2 R /d6\2 , . d9 ,~ 9 , o 

w - R(-) - {v + 2*)R-- { p + a 2 )R- 9l R* 

= [V lC os 2 (k£) + f 2 cos 2 (7fc£ + </>)]#, (4) 
d 2 9 dRdO , ,dR 
V + 2 ^ + (U + 2a) ^°- (5) 
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In the formulas we have used the dimensionless vari- 
ables and parameters £ — fc C, v — 2mvi,/hko, (3 = 
hp/E r0 , a = a/ko, Vi = Vi/E r0 , V 2 = V 2 /E r0l E r0 = 
h 2 k 2 /(2m) and gi = Aa s /(kol 2 ) with k being the unit 
of wave vector and l r = y/K/ '(raw,) the radical length of 
harmonic oscillator. Integrating Eq.© yields dd/d£ — 
C/R 2 — (v/2 + a), where C is the integration constant. 
Applying this to Eq. ([5]) , we arrive at the decoupled equa- 
tion 
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va — p. It is very hard to find the 
exact solution of this equation. However, when the 
driving strengths are weak enough, we can treat the 
chaotic system by the direct perturbation approach [26j 
and Melnikov- function method [2 71]. 

For simplicity we consider the C = case which leads 
Eq. (6) to the driven Duffing equation, whose chaotic 
features have been extensively studied for the single lat- 
tice case [10, [2l|. In the double lattice case, we expand R 
to the first order, 



^(0 = ^(0 + ^1(0, l#o|>|#i| 
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and insert it into Eq.(6) with C = 0, obtaining the 
zeroth-order and first-order equations 



£ + DRo- g x Rl = 0, 
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e(0 = [V 1 cos 2 (kO + V 2 cos 2 (jk^ + ^)}Ro. 



If the atom-atom interactions are attractive, the system 
has a negative s-wave scattering length to make g\ < 
such that Eq.® with a negative D has the well known 
homoclinic solution 
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Here £o = ^0(^0 + ^l^o) is the combination of the initial 
time to and boundary coordinate xq, cq denotes an inte- 
gration constant determined by the boundary and initial 
conditions. Obviously, the zeroth-order number density 
i?o(0 is just a bright soliton solution. Applying Eq. fTT)]) 
to Eq. © , we construct the general solution of Eq. © in 
the integral form[2t| 

Ri = h t f e(0 <%-f [ he(£) d£, (11) 

JP JQ 



where P and Q are the integration constants, / = dRo /d£ 
and h = / J .f~ 2 d£ are two linearly independent solutions 
of Eq.© for e(£) = 0, 
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x {eV 17 ^^ + c ) - 4coth[V^(£ + c )] 
+ sinh[2V^D(£ + c )]}. 
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Generally, the corrected solution (fTTj) is unbounded, be- 
cause of the unboundedness of h at £ = 00. However, 
using the l'Hopital rule, we can easily obtain the suffi- 
cient and necessary boundedness condition [2t| 



7±= lim / /e(0de = 0. 
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It is worth noting that the relation between the Mel- 
nikov function M(cq) and Eq. (fT3|) is M(co) = /+—/- = 
/-» / d ^ Combining e(£) in Eq.® with R in 
Eq. (fT0|) and / in Eq. (fT2|) to the integrand, we yield the 
well-known Melnikov chaos criterion [23 



M(c ) = 
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sin(2co7/c - 2(j>) = 0, (14) 



which indicates the existence of chaos for some Co values. 
Under the conditions of Eqs. fTS"]) and (jTJJ) we can call 
Eq. pTj) the "chaotic solution". Thus inserting Eqs. (10) 
and (11) into Eq. (7) produces the chaotic bright soliton 
solution which is the superposition between the soliton 
and chaotic states and propagates with the velocity of 
traveling superlattice. 

Chaos Probabilities and Regions The Melnikov func- 
tion is a periodic function of Co for the fixed parame- 
ters and rational number 7, so only the discrete zero 
points Co = Cot for i = 1, 2, • • • satisfy the chaos crite- 
rion M(coi) = 0. Nevertheless, Co is an integration con- 
stant depending on the initial and boundary conditions 
and cannot be accurately set in experiment. When the 
stochastic initial and boundary conditions are applied to 
the numerical or experimental investigations, from Eq. 
(10) we know that Co takes coi values with only a certain 
probability which is just the probability from soliton to 
chaos. Defining this probability as the chaos probability 
P, it is clear that the P value should be proportional to 
the number n of coi in one period of M(co) and is always 
smaller than 1. Taking the chaos probability of single 
lattice case as the referential one Pq, we have the rela- 
tion P = hPq/2 between P and n, since n is equal to 2 
in the single lattice case as in Eq. (14) with V 2 = 0. 

Now let us see how the chaos probability depends on 
the parameter regions. As a simple example we firstly 
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zero point and the number n = 2 of zero points of M(cq) 
agrees with that of X\(cq). 



FIG. 1: Plot of the chaotic regions of V 2 versus k for the dimcnsionlcss 
parameters 7 — 2. D — — 2 and Vi — 0.04. 



consider the case 4> = and 7 = N > 2 with N integer. 
We rewrite Eq. (fl4| in the form 



M(c ) = rjX 1 (c )X N (co)=0, 



(15) 



where rj = 2k 2 ir/gi, Xi(co) = sin(2cofc), and 
X N {c ) = 7iCSch(fc7r/v /= S) + F N (c ) with F N (c ) = 
y 2 JV 2 csch(^) Sgggl . Clearly, the X^cq) is a 
periodic function of c but with different number 
of zero points in different parameter regions, which 
can lead to different number n of coi and different 
chaos probabilities. This can be easily evidenced in 
the case N = 2. The super lattice of Eq.(l) with 
7 = 2 has be widely studied 24]. In such a case, 
the Xn(co) becomes X 2 (cq) = ViCscli(A;7r/ \J— D) + 
8V^csch(2fc7r/\/— D) cos(2cofc) and M(cq) has period tt. 
Obviously, two parameter regions with different zero 
points of X 2 (co) exist and the boundary between them is 
given by 



V 2 = Vi 



2h 



Vicsch(fe7r/-\ 
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As an instance from Eq. (16) with D = —2, Vi = 0.04 the 
boundary curve on the parameter plane (k, V2) is plotted 
as in Fig. 1. The region A above the boundary curve 
is associated with V 2 > V 2 b and the region B below the 
boundary curve corresponds to V 2 < Vib- 

In order to show the different numbers n of zero-point 
CQi in one period of M(cq) for the different parameter 
regions, we make the sketch maps showing the zero points 
of functions Xi(co), X 2 (c ) and M(co) as in Fig. 2. For 
the single lattice case with V 2 = from Fig. 2(a) we see 
that X 2 (cq) = constant, M(cq) and Ai(c ) have n = 2 
same zero points, indicating the chaos probability P — 
nPo/2 = Pq. In Fig. 2(b) with parameters of region 
A, we exhibit that both Xi(co) and X 2 {cq) have two 
zero points in range Co € [0, tt) respectively, and all the 
zero points are not coincident. This means n = 4 and 
the chaos probability P — nPo/2 = 2Pq in region A. 
On the boundary between regions A and B, from Fig. 
2(c) we can see that two zero points coincide among the 
three zero points of Xi(cq) and X 2 (cq) in one period. 
Therefore, we have n = 2 and P = Pq on the boundary 
curve. The same chaos probability appears in region B, 
which is illustrated by Fig. 2(d), where X 2 (cq) has no 
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FIG. 2: Sketch maps of functions X\{cq) (dashed curve), Xz^cq) 
(doted curve) and M(cq) (solid curve) versus Co for the parameters 
7 = 2, D - -2, gi = -0.5, k = 1, Vi = 0.04 and (a) V 2 = 0, 
(b) V 2 = 0.08 > V 2b , (c) V 2 = 0.01cosh(7r/ v / 2) = 0.0466 ss V 2b , (d) 
V 2 — 0.01 < V2b, where for each different function the different ampli- 
tude scale is adopted for showing the zero points. 

The above results can be classified as the two cases: 
Case 1. In the double chaotic region A with pa- 
rameters obeying V 2 > V 2 b, the chaos probability reads 
P = 2P ; 

Case 2. In the chaotic region B and on the bound- 
ary curve with parameters obeying V2 < V 2 b, the chaos 
probability equates the referential one, P = Pq. 

To numerically confirm the analytical results, we use 
the MATHEMATICA code 

T = 7r, e[{Rnewvnewj] := {R[T], v[T]} / .Flatten 
lNDSolve[{R'lZ} == == <?ii?[£] 3 - DR[t] 

+ {V 1 Cos[k£] 2 + V 2 Cos[ykt] 2 )R[£\,R[0] == Rnew, 
v[0] == vnew}, {R, v}, {^, 0, T}}]; Do[pici = ListPlot 
[Drop[Nestlist[e, {Random[Real, {-2.1, 2.1}], 
Random[Real, {-2.2,2.2}]}, 3000], 100], {i, 1, 500}] 

to make two groups of Poincare sections on the equiva- 
lent phase space (i?, R^) for the parameters used in Fig. 
2(b) and Fig. 2(d) and the random initial conditions 
{R(Zo) G [-2.1,2.1], f^(fo) G [-2.2,2.2]} associated 
with a suitable range of Co- Each of the groups contains 
five hundreds of Poincare sections. For the parameters in 
region A, 180 chaotic attractors are found, and the num- 
ber of chaotic attractors is 85 for the parameters in region 
B. The numerical results show that in the chaotic region 
B the chaos probability reads P = P = 85/500 = 0.17, 
and in the double chaotic region A the chaos probability 
becomes P = 180/500 = 0.36 w 2P . They well agree 
with the analytical assertion. The small difference exists 
between the analytical and numerical results, because of 
that the strict chaos probability requires more samplings 
of the Poincare sections than the used ones. 

Then, we take 4> — and 7 = = 3, 4, 5 into account 
and find that for any the parameter space can be di- 
vided into several parts with different chaos probabilities. 
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The lowest chaos probability is Pq and the highest one 
arrives at NP . Moreover, the number of chaotic regions 
increases with N. 

Control of Chaos Probabilities Controlling the chaos 
probability to zero is important for some cases one re- 
quires suppressing chaos. On the other hand, for the 
secure chaos-based cryptosystems[22|, HH one calls for 
higher chaotic probability. Given the above-mentioned 
results, we can decrease or raise the chaos probability by 
adjusting the parameters into different chaotic regions. 
In order to suppress chaos, at first we let the parameters 
enter the region of the lowest chaos probability (e.g. the 
above region B), where only Xi(co) = sin(2coA;) has zero 
points. For such a parameter region if a chaotic state 
is observed experimentally, which implies 2cofc = i~n for 
% = 1,2, ••• and the first term of Eq. (14) vanishing. 
Then we can adopt the two different methods to elimi- 
nate chaos: 

Method 1. We take = and vary 7 from an integer 
to an irrational number suddenly such that the second 
term in Eq. (14), sin27Cofc = sin 1777, does not vanish 
and the Mclnikov chaos criterion M(co) = cannot be 
fulfilled, namely n = and P = nPo/2 = 0. 

Method 2. Let 7 be equal to 2 and suddenly change 
4> from zero to nonzero, say 7r/4, the nonzero Melnikov 
function can also be obtained, that leads the chaos prob- 
ability to zero. 

The operation of method 1 can be performed by switch- 



ing off the second lattice with 7 = integer suddenly and 
switching on the third lattice with 7 = irrational number 
simultaneously. Similarly, we can complete the opera- 
tion of method 2 by a new lattice of phase difference 7r/4 
with the second lattice. The analysis offers an effective 
technique for suppressing spatiotemporal chaos experi- 
mentally 

In conclusion, we have investigated a BEC system 
loaded into a weak moving optical superlattice created by 
interference between two lattices of different depths and 
wave vectors. The superlattice separates the chaotic re- 
gion into several parts with different chaos probabilities, 
which are studied analytically and numerically. More- 
over, for a fixed first lattice, the modulation of the sec- 
ondary lattice can transform the chaos probability to zero 
or higher one. Our results suggest a feasible method 
for suppressing or strengthening chaos experimentally. 
When the zero traveling velocity vl = is set, the re- 
sults fit the corresponding static BEC system. The used 
method could also be easily extended to the controls of 
other chaotic systems. 
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